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ABSTRACT 



Simultaneous measurements of wave height and two orthog- 
onal water particle velocities were made at six elevations 
in nineteen meters of water using a penetrating v;ave staff 
and an electromagnetic flowmeter. Moderate sv/ell and low 
wind conditions prevailed during the experiment. The mea- 
sured wave-induced velocities were tvjo to four percent 
greater than those predicted using linear xvave theory. 
Coherence of the wave height and wave-induced velocities in 
the significant energy-density range was computed to be over 
0.85, Indicating that the motion was almost totally wave 
induced. At higher frequencies it v^as apparent that the 
motion was primarily turbulence. Phase spectra computed for 
the measured v;ave heights and orbital velocities compared 
very well with linear theory. Measured frequency distribu- 
tions were compared to both Gaussian and Gram-Charlier 
distributions using the chi-squared goodness-of-f it test. 
Qualitatively, the Gram-Charlier distribution gave the 
better fit to the' data. 
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I . INTRODUCTION 



A. BACKGROUND 

Measuring Instantaneous water particle velocities in 
the presence of wave-induced motion has been difficult due, 
primarily, to a lack of instrumentation. This is particu- 
larly true for measurements made in the field. Until re- 
cently, there has been no adequate instrument available to 
measure rapidly fluctuating flows found beneath wind v;aves 
and swell. As vrell as the usual requirements of reliability 
and accuracy, such an instrument must also be rugged, have a 
low response time, and have good directional measuring 
characteristics . 

Earlier attempts at correlating v/ater height and particle 
motion have generally been accomplished in shallow water next 
to the shore, off piers, or in estuaries. Historically, 

Inman and Nasu (1956) made one of the earliest attempts at 
field measurements of wave-induced viater particle motion, 
using a force meter in shallow v/ater and compared their 
results to solitary wave theory. Miller and Zeigler (1964) 
used both an acoustic meter (based on the doppler shift 
principle) and an electromagnetic flowmeter to measure water 
particle motion in the surf zone. They compared their 
measurements to higher order wave theory and obtained 
qualitative agreement. 

Nagata (1963, 1964a, 1964b) used an electromagnetic 
flowmeter in a series of experiments to measure water 
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particle velocities in the surf zone and to determine wave 
directional spectra. Bov.’den and V/hite (1966) and Simpson 
(1968) performed field measurements using an electromagnetic 
flowmeter developed at the National Institute of Oceanography. 
Like Nagata, Bov/den and White (1966) used water particle 
velocities to determine the directional spectrum of waves. 

Simpson's experiment vjas performed at an open coastal 
site in vjater approximately 6 meters deep. The v/aves, 
primarily v;ind generated sea waves, v;ere developed over a 
limiting fetch of 210 km. Simpson employed spectral methods 
in the analysis of the vjaves and of the v:ave-induced water 
particle velocities. His results were consistent with 
linear wave theory. 

Seitz (1971) carried out a similar experiment in that he 
measured sea surface displacement and vjater velocity compo- 
nents in an estuary at a depth of one meter. He compared 
the measured velocity spectra to a derived velocity spectra 
obtained from direct measurement of the surface wave field 
using linear theory. Seitz had to contend with a shift in 
the observed frequencies, the magnitude of v;hich was a 
function of the tide induced velocity. Since at the fre- 
quency of peak amplitude the error introduced was approxi- 
mately 6.7^, he did not correct the data. The total v/ater 
particle motion at depth can be described as the sum of the 
turbulent velocities and the surface wave induced motion. 

Seitz attempted to separate the two based on the slope of 
the spectrum in the Inertial sub-range. 
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The lack of Instrumentation has not been the only 
hindrance to wave height and velocity studies. Most research 
has been conducted at shallov; v/ater inshore sites because of 
the general unavailability of a deep water, exposed instru- 
ment platform. Shonting (1966) used the Buzzard's Bay 
Entrance Light Station on the South Coast of Massachusetts 
to obtain a twenty meter vmter depth at an open ocean site. 

He used ducted current meters to measure the orthogonal 
components of the water velocity under waves. His experiment 
was designed primarily to evaluate these meters. An inherent 
difficulty with any type of ducted meter is their poor 
directional measuring characteristics when the angle between 
meter axis and the flov; is large. 

Bordy (1971) made measurements in 20 meters of v;ater 
from the Naval Undersea Research and Development Center (NUC) 
Oceanographic Research Tower located approximately one mile 
off Mission Beach, San Diego. He compared v/ave height and 
water particle velocity using spectral analysis, and found 
a general agreement with linear vjave theory except for 
measurements of phase spectra at the deepest measurements of 
ten meters. The phase shift encountered v^as unexplained in 
the phase angle error analysis. During the experiment, the 
NUC tower's permanently installed wave staff system was used 
for water level sensing. The flowmeter was removed horizon- 
tally from the v;ave staff by approximately five feet. It is 
felt that due to this displacement, an error was introduced. 
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During 8 June 1972 records of orthogonal v/ater particle 
velocity and instantaneous water elevation were again 
measured at the NUC tov;er. The data vjere taken in conjunc- 
tion v;ith an experiment involving the investigation of vave 
induced fluctuations in acoustic phase shifts, temperature, 
and salinity. 

The results of the wave and water particle motion 
measurements are described here. An electromagnetic flow- 
meter was used to gather the wave particle velocity data 
at several depths. These measured values of v;ave-lnduced 
motion were compared v;ith values predicted using linear 
theory and the wave height record. 

B. RESEARCH OBJECTIVE 

The purpose of this investigation v;as four fold: 

1. To measure simultaneously wave height and two 
orthogonal wave-induced v;ater particle velocities. 

2. To perform a spectral analysis on the data to obtain 
its statistical properties. 

3. To spectrally compare the data to first order wave 
theory. 

4. To determine the statistical geometry of both the 
sea surface and the measured velocities. 
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II. EXPERIMENTAL SITE AND INSTRUMENTATION 



A. EXPERIMENTAL SITE 

The experiment was conducted at the Naval Undersea 

Research and Development Center Oceanographic Research Tov;er 
\ 

at San Diego, California, 8 and 9 June 1972. The NUC tov/er 
is located approximately one m.lle off Mission Beach, San 
Diego in a mean water depth of 19 meters. It is constructed 
of steel and concrete and is permanently imbedded in the sea 
floor. Figure (1). The forty foot square tower has three 
levels. The first is used primarily for gaining access, 
while the second houses machinery and provides work space 
for placing instruments in the water. The thir'd, which is 
enclosed, is used for the housing of electrical equipment 
and also provides sleeping space for six men. The legs of 
the tower slope five degrees. A pair of tracks is mounted 
on the north, south, and v;est sides for lov;ering scientific 
instrument packages into the sea to any depth above the 
ocean floor. The v;estern side of the tower was used for the 
experiment because of its exposure to the predominant swell. 

B. INSTRUMENTATION 

Water particle velocities v;ere measured with an Engineer- 
ing Physics Company v;ater current meter Model EMCM-3B. The 
sea surface elevation v;as measured with a Baylor Company 
wave staff system Model 13528R. The penetrating wave mea- 
suring system gives a direct measure of the sea surface 
elevation. 
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Figure 1. Naval Undersea Research and Development Center 
Oceanographic Research Tower, San Diego, California 
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The model 13528R Baylor Company VJave Staff System is an 
Instantaneous v;ater level measuring device consisting of 
two tensioned 1/2", 6x19 IWRC 30^ SS wire ropes. Together 
with these lengths of stainless steel wire rope, the trans- 
ducer produces an electrically linear direct current output 
that is proportional to the amount of v;ave staff above a 
short circuit produced by the water surface. System accuracy 
is stated as within one percent of actual height and the time 
constant is less than 0.06 seconds. 

As received from the manufacturers, the v;orking length 
of the staff was 50 feetj however, in order to take water 
particle velocity measurements near the surface directly 
under the staff, it was shortened to 20 feet. After 
shortening, the staff v.’as recalibrated. The resulting 
calibration curve is shown in Figure (2). The rails on 
which the Instrument cart rides slopes with the legs of the 
tower. This necessitated the rigging of an aluminum "I" 
beam with a trolley to support the Baylor wave staff system 
directly over the electromagnetic flov/meter. The wave staff 
was tensioned by allowing the weight of the instrument cart 
and frame to be exerted on the lower end of the wave staff 
cables. The estimated tension was approximately 800 pounds. 
The instrument and installation are shovm in Figure (3). 

The EMCM-3B electromagnetic current measuring device 
measures two orthogonal components of water particle velocity 
through a range of zero to five M/SEC. Other characteristics 
of the flow meter are: 
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Figure 2. Calibration Curve for the Baylor Company Model 13528R V/ave Staff System 




Figure 3. Model 13528R Baylor Company Wave Staff System 
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Output Voltage: ±5 volts full scale 

Electrical Time Constant: 0.2 seconds or longer, 

adjustable by user 

Error Band: the largest of 

1) ±1 percent of full scale 

2) ±5 MM/SEC 

3) 5/(T^) MM/SEC 

where T is the time constant 

Transducer Size: 2 3 /^" O.D. x 36" L 

The operation of the unit is based on Faraday's principle of 
magnetic induction. A current is induced when a conducting 
fluid flows around the probe through an alternating magnetic 
field generated by a toroidal magnet inside the instrument. 
The equation describing the production of current is: 

g = u X g, 

where 2 is the induced voltage, u the velocity of the 
conducting medium, and § the magnetic field. The Induced 
electromotive force is linearly proportional to the amount 
of fluid floviing through the magnetic field. The magnetic 
field, B, is most intense in the near vicinity of the probe 
and decreases according to the inverse square of the distance 
from the instrument. Accordingly, the flow velocity near 
the probe is weighted most heavily in its contribution to the 
measured velocity. Qualitatively, only that water which is 
within a distance of tv;o to three probe radii contributes 
significantly to the flow induced voltage. 
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The flow meter was originally calibrated under steady 
flov; conditions by the manufacturer. Different response 
characteristics are to be expected under unsteady flow 
conditions; consequently, the probe was recalibrated by 
oscillating it in a viater tank. The equipment schematic is 
shown in Figure (4). The carriage on which the flov.' meter 
probe was mounted travels back and forth on rails. It is 
driven by a variable speed electrical motor. The peak 
carriage velocity was calculated from the tangential 
velocity of the motor arm. The ratio of carriage velocity, 
V, , to the velocity measured by the instrument, V , v/as 
found for different angular velocities. This was done for 
each pair of electrodes by orienting them parallel to the 
flow in turn. 

The instrument's response was vei’y linear with respect 
to frequency up to at least 0.5 Hz. Beyond this frequency 
the cart upon which the probe v;as mounted became unstable. 
Thus 0.5 Hz was the maximum frequency Investigated during 
calibration. This limitation v;as not deemed a debilitating 
factor since 0.5 Hz was above the upper frequency limit of 
significant wave energy. The results of the calibration are 
shown in Figure (5) and the Instrument and installation in 
Figure (6). 
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Figure 4. Schematic of Calibration Equipment for Flowmeter 
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Figure 5. Calibration Curve for the Engineering Physics 
Company Model EMCM-3B Flowmeter 
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Figure 6. Engineering Physics Company Model EMCM-3B 
Electromagnetic Flowmeter Installed in Position (arrows) 
on the Instrument cart 
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III. EXPERIMENTAL PROCEDURE 



Sea surface elevation and tv/o orthogonal components of 
water particle velocity were measured simultaneously. The 
cart on v/hich the flowmeter was attached v;as raised and 
lov/erod so that measurements could be made at var j ous depths. 

Data from the flovmeter and v;ave gauge v;ere recorded on 
a Sangarno Model 3500 fourteen channel I'M tape recorder 
during the ten data collection periods. The data chosen for 
analysis v/as obtained during the first six recording periods 
on 8 June 1972. Each processed run was about 20 minutes 
duration. After initial processing, the other four data 
collection periods v;ere deemed too short to be statistically 
representative . 

Summarized in Table I below are the velocity components, 
water depth, depth of the instrument, and record length for 
each of the processed runs. Instrument depth refers to the 
depth of the flov;meter below mean sea level. The u-component 
of velocity was oriented north-south approximately parallel 
to the shore, vjhile the v-component was measured east-v;est. 
Correct flovmeter orientation V7as obtained using a plumb-bob 
and level. The v;ave gauge was always maintained directly 
over the flowmeter. As the cart was lowered the top of the 
wave staff was pushed further avzay from the tower to compen- 
sate for the five degree tilt of the tracks to which the 



cart was fastened. 
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The v/ave staff’s perpendicularity to the horizontal v;as also 
ascertained by the use of a level. The v/ave staff installa- 
tion is pictured in Figures (7) and (8). 



TABLE I 



RUN 


DEPTH 
(meters ) 


DEPTH OP 
INSTRUMENT 
(meters ) 


COMPONENT 

OP 

VELOCITY 




LENGTH OF 
RECORDS 
(minutes ) 


1 


19.15 


7.04 


u & 


V 




24. 


96 


2 


19.22 


14.62 


u & 


w 




22. 


61 


3 


19.25 


12.64 


u & 


v; 




21. 


33 




19.28 


10.69 


u 8c 


V7 




24. 


85 


5 


19.32 


8.66 


U 8: 


W 




18 . 


88 


6 


19.36 


5.70 


u 8c 


V7 




19 . 


52 


The sea surface was in general disturbed 


by 


swell coming 


from 


the v;est- 


southwest. The sky was 


overcast 


with 


a low 


stratus layer 


and little wind. 


less than 10 ; 


knots generally. 


Sea 


surface temperature was obtained using a 


bucket 


thermom- 


eter 


. Table II summarizes the 


environmental 


data. 
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TIME OP 


AIR SEA 




WIND 




CLOUD 


SV/ELL 


RUN 


START OF 


TEMP SURF 


DIR 


VEL 


COVER 


. DIR 




RUN (PDT) 


(°C) TEMP i°C) 


(°T) 


(knots ) 


(%) 




1 


lil05 


19.7 19.20 


240 


8 




100 


WSW 


2 


1600 


19.6 18.95 


250 


10 




100 


WSW 


3 


165 ^ 


19.2 18.20 


240 


9 




100 


WSW 




1800 


18.7 19.10 


225 


6 




100 


WSW 


5 


1830 


18.2 18.99 


225 


7 




100 


WSW 


6 


1922 


17.9 18.77 


235 


4 




100 


WSW 
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Figure 7. Model 13528R Baylor Wave Staff Installation 
Showing "I" Beam and Trolley for the Top of the System. 




Figure 8. Pulley Attachment Point for the Bottom of the 
Baylor Wave Gauge (arrow). 
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Table III summarizes the temperature profile data taken 
at the beginning of each data collection period. This data, 
which shov/s a remarkable variation in layer depth and 
Intensity of temperature gradient, v;as obtained using a 
mechanical BT. The NUC tov;er v;as originally constructed to 
study Internal waves which are prevalent in the area. The 
variability in temperature profiles is attributed to Internal 
waves. A typical temperature profile is shown in Figure ( 9 ). 



TABLE III 



RUN 


BT 
SURF 
TEMP ( ° C ) 


TOP OP 

ist layer 

(°C/m) 


BOTTOM OP 

ist layer 

("C/m) 


TOP OP 
2 i^‘d LAYER 
(°C/m) 


BOTTOM 

TEMP 

(°C) 


1 


22.8 


22.2/7.31 


16.6/12.19 


xxxxxxxxxx 


15.55 


2 


22.5 


22.2/3. OM 


19.9/12.80 


17.2/15.85 


16.66 


3 


23.3 


23.3/6.08 


21.6/9.1^ 


20 . 5 Zl 6 .i 46 


18.33 


ij 


21.1 


21.1/^. 87 


19 .^/ 6.71 


18.3/17.68 


16.66 


5 


22.2 


21. 6/5 . ^9 


19.9/10.67 


17.7/17.06 


18.33 


6 


19.9 


19.9/5.79 


16. 1 / 15 . 2^4 


xxxxxxxxxx 
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The analog PM recordings were digitized using a hybrid 
system consisting of a Scientific Data Systems (SDS) Model 
5000 analog computer and a Xerox Data Systems Model 9300 
digital computer. For comparison purposes the data were 
transcribed onto a rectilinear eight channel Clevite strip 
chart recorder after a tenfold amplification through the 
(SDS) 5000 analog computer. An exam.ple of the transcribed 
data, after digitization, is shown in Figure ( 10 ). After 
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Figure 9 . Typical Temperature-Depth Profile, 
Run Number 4 . 8 June 1972. 
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Figure 10. Transcribed Analog Data after Digit Izln 
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digitizing onto the seven track tape used by the 9300 system, 
the data v;ere transcribed onto a nine track tape using the 
Naval Postgraduate School’s IBM Model 360 digital computer. 
This computer was also used for the statistical analysis. 
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IV. TIIEORY 



A. LINEAR WAVE THEORY 

The elevation of the surface n(t) can be described as 
the superposition of an infinite number of sinusoids of the 
form: 



n(t) = I a cos (i^ -x-a- t+e ) = I n 
n=l ^ ^ n=l ^ 



(^. 1 ) 



v;here x is the horizontal Cartesian coordinate, t is the 
time, is a horizontal vector wave number, the frequency 
related in linear theory to k by : 



a = Kk tanh k h 
n ^ n n 



(^. 2 ) 



where g is the acceleration of gravity, a^ is the amplitude 
and the phase angle. ri(t) represents the sum total of 
all component vmvelets. Summing in the manner of Equation 
(^.1) implies a linear system and restricts the analysis to 
the use of linear or Airy wave theory to describe the vjave 
induced motion. 

Linear v/ave theory can be used to predict the wave 
induced water particle velocities. The equations for the 
horizontal and vertical velocities respectively are: 



u 



~ a^gk^ cosh k (h+z) 

(t) = I cos(K„ *x-a„t+e„) 



n=l 



a sinh k h 
n n 



n 



n 






or 



» gk cosh k (h+z) 
__ r n n 

" a sinh k h 

n=l n n 



(^. 3 ) 



n 
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and 

w(t ) 



00 



1 

n=l 



slnh k^(h+z) 



a sinh k h 
n n 



sin(K • x-a t+e ) 
n n 



00 



= I 

n=l 



gk^ slnh k^(h+z) 

^ sinh k h 
n n 




e=--90° 






where h is the total v;ater depth and z is the depth of 
interest measured positively upward from the still v;ater 
level. 

The terms 




cosh k^(h+z) 
sinh k^h 



and 

gk^ sinh k^(h+z) 

d sTnh k h 
n n 






(4.6) 



represent the transfer functions relating the horizontal and 
vertical velocities to n(t). 

B. LIMITATIONS OF FIRST ORDER WAVE THEORY 

In the formulation of linear wave theory, the boundary 
conditions are linearized in order to obtain an analytical 
solution. In the linearization, it is assumed that the 
amplitude, a, is small compared to the wave length, L, that 
is, a/L << 1. 

Higher order solutions to the boundary value problem, 
generally obtained by perturbation analysis, give a better 
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representation of a constant profile wave. Hov;ever, the 
nonlinearities introduced in the solution precludes their 
use where the principle of superposition is invoked. This 
in turn eliminates their generalization to a stochastic 
process . 

Dean (1968) analyzed the various v;ave theories with 
respect to hovr well they fitted both the dynamic and kine- 
matic free surface boundary conditions. He showed that 
first order theory provided a good solution to the v/ave 
equation for vraves of permianent form having small steepness 
in intermediate water depths (the assumption of linear 
theory stated above). 

The conditions prevailing during this experiment fell 
v;ithin these limitations. The rms v/ave height averaged 0.27 
meters with a mean period of I6 seconds corresponding to a 
deepv;ater wavelength of 209 meters in 19 meters of water. 
Hence, if linear theory provides a reasonable approximation 
for the average wave conditions, the superposition of a 
large number of components, each described by linear theory, 
should provide an even better description. 
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V. THEORETICAL CONSIDERATIONS OF DATA PROCESSING 



BY STATISTICAL ANALYSIS 



The digitized records v;ere spectrally analyzed using 
the methods of Blackman and Tukey (1958) and utilizing the 
IBM Model 360/67 computer at the Naval Postgraduate School. 

A. COMPUTATION OP THE ENERGY-DENSITY SPECTRUM 

The covariance function ({) ,(x) of time series random 
data describes the general dependence of the values of the 
data at one time on the values at another time for the same 
data set. In integral form 

T 

i / x(t) x(t+T) dt (5.1) 

vjhich is the average of the product of the value of a data 
point x(t) with that of an adjacent point at time t later. 

The resulting average product will approach the exact auto- 
correlation function as T approaches infinity. 

VJhen using a covariance function generated v;ith a finite 
data set, it is necessary to apply a lag window to <P ^ (x). 

A X 

A Parzen vrindovj, as used in this analysis, has the desirable 
property of no negative side lobes thus eliminating any 
possible numerical instability v;hich might arise due to 
resonance effects of the side lobes. This is particularly 
critical in computing cross-spectra. The Parzen window is 
defined as: 



3 ^ 



P(x) = 1 - 6 (t/T^)^ + 6(t/Tj^)^ t = 0,1,2, ... ,Tj^/2 

= (1 - (t/T ))^ X - T /2 + 1,...,T 

m IB m 

= 0 X = T 

m 

(5.2) 

vihere T Is the total lag time. 

The energy-density spectrum as computed from the Fourier 
transform of the covariance function viith the Parzen viindow 
applied to correct for finite record length is 

oa 

= / P(x) dx , (5.3) 

CO 

where f denotes a particular frequency. 

The energy-density spectrum is the statistical estimate of 
the energy-density distribution in the various frequency 
bands comprising the spectrum and the area under an energy 
density spectrum is equal to twice the variance of the random 
variable. With respect to the wave height spectrum>the area 
is proportional to the total potential energy-density of the 
sea surface elevation, vihile the area under the velocity 
spectrum is proportional to the total kinetic energy-density. 

B. COMPUTATION OP THE CROSS-SPECTRAL DENSITY 

The cross-covariance function is similar to the covariance 
fucntion with the exception that it compares two time records 
lagged by a time interval x. Thus, in integral form 
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<j)j.y(T) = i / x(t) y(t+T) dt . 

If X and y are statistically Independent and the mean of 

both X and y is zero, <p (t) = 0 for all f. If not, <{> (x) 

xy xy 

is equal to the product of the means. The lag window 
considerations for a finite record length are identical to 
those for the covariance function. The cross-spectral 
density is therefore the Fourier transform of the Parzen 
window lagged cross-covariance function 

oo 

$^y(f) = / P(t) dx ( 5 - 5 ) 

The cross-spectrum can be defined in terms of its real 
and imaginary parts since, unlike the energy-density spectrum., 
it is not an even function. Thus 



$ (f) = C (f) - iQ (f) 

xy^ ^ xy^ ^ ^xy^ 



(5.6) 



where the real part of $^y(f) is the co-spectrum 



C^y(f) = 2 / [<|).^^(x) + cti^^(-x)] cos(2TTfx) dx , (5-7) 



xy 



xy 



and the imaginary part the quadrature spectrum 



Q (f) = 2 / [<|) (x) - (j) (-x)] sin(2irfx) dx . (5.8) 

xy Q xy Ay 
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C. COMPUTATION OP COHERENCY 



A useful quantity arising out of the computation of the 
cross- and energy-density spectra is the coherency. 
Coherency is defined as 



(f)]‘ 

xy 



- - 

^ (f) (f) 



XX 



yy 



(5.9) 



where 

0 1 Y^xy^^^ - ^ * 



( 5 . 10 ) 



2 

V/hen Y v^ 7 (^) “ records x and y are 

xy 

fully coherent and statistical dependence can be assumed. 

2 

On the other hand, if y xy^^^ ~ ^ ^ 

statistically Independent. Bendat and Plersol (1966) 
attribute coherency values less than unity to three possible 
causes: the system is not linear, extraneous noise is 

present in the measurements, or the output response is due 
to more than one input function. 

D. COMPUTATION OF PHASE ANGLE 

The cross-spectral phase angle represents the average 
angular difference by which the cross-correlated components 
of y(t) lead those of x(t) in each spectral frequency band. 
The phase angle can be computed using the co- and quadrature- 
spectra 

e^^y(f) = arctan [Q^y(f) / C^y(f)] . 



(5.11) 



E. PROBABILITY DENSITY FUNCTION OF THE SEA SURFACE 
AND INDUCED V/ATER PARTICLE FLUCTUATIONS 

The probability density function, p(x), of random data 

describes the probability that the data will assume a value 

within some defined range at any instant of time. 

Mathematically 

, V _ lim Prob[x<x(t )< (x+dx) ] _ lim _1_ rlim 

dx->0 dx dx->-0 dx T 



where T^ is the amount of time that x(t) falls in the 
range (x, x+dx) during an observation time T. 

For this experiment, the sea surface displacement u’as 
assumed to be the sum of an infinite number of independent 
components of random phase. Under these conditions the 
measured frequency distribution of the sea surface displace- 
ment tends to be Gaussian. For zero mean the Gaussian 
probability density function is given by 



PqCC) 



i- 

2it 



(5.13) 



where C is equal to the value of the measured quantity 
divided by the standard deviation. However, if the non- 
linearities of the boundary condition become important, 
departures from Gaussian are expected. Longuet-Higgins 
( 1963 ) showed for "moderately" non-linear swell in deep 
water that the probability distribution of n is better 
described by the Gram-Charlier probability density function. 
Pg^(^) is given by: 
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= Pq(5) [1.0 



^ + -j7^ Hg) + ...]. (5.1^) 



Here the ' s are cumulants with ra^ being the skev/ness of 
the data record, the kurtosis minus three, and the 
Hermit e polynomials of degree n. The Hermite polynomials 
are generated from 



TT _ n(n-l) ^ . n(n-l) (n-2) (n-3) 

n “ 1 ! 2 2 ! , ^2 

Therefore 

H 3 = - 3 , 

2 i . 2 

H|^ = - 6 ^^ + 3 , 

and h o 

Hg = - 155^ + ^5C - 15 . 



(5.15) 



( 5 . 16 ) 



If the distribution of 5 is Gaussian, the skewness is zero 
and the kurtosis is equal to three v;hich results in these 
and all higher order cumulants being zero making 

Kinsman (i 960 ) made measurements of wind generated v;aves 
and showed their distribution was only approximately Gaussian 
and better approximated by the Gram-Charlier distribution. 

He found the distribution exhibited a positive skewness and 
attributed this skeviness to the physical fact that waves have 
a distinct tendency to form short sharp crests and long 
shallov; troughs. 



39 



i 



F. GOODNESS-OF-FIT TEST 



The chl~square' goodness-of-flt test is used to compare 
a measured probability density function (PDF) to a theoret- 
ical one. A statistic v/ith an approximate chi-squared 
distribution is calculated as a measure of the discrepancy 
betvjeen an observed PDF and the theoretical PDF. The 
chi-squared parameter is: 



where f^^ is the observed frequency of occurence, F^^ the 
theoretical frequency of occurence, and J is the number of 
intervals into which the frequency distribution is divided. 

If a calculated distribution miatches exactly a theoretical 
distribution the chi-squared parameter is zero. 

Two methods can be used to calculate the chi-squared 
parameter for the goodness-of-flt test. In the first method, 
the ^-axis of the PDF is divided into J Intervals of equal 
size. The f ' s in this case represent the number of data 
points in each Interval betv/een divisions and f^ varies 
across the distribution. In the second method, the PDF is 
divided into J intervals of equal probability density, each 
containing approximately an equal number of data points in f^. 
V/illlams (1950) showed that to give the test sufficient 
power the number of Intervals, J, should be given by 




(5.17) 



J = 2 




( 5 . 18 ) 



^10 



v/here n is the number of time samples. C, under the null 
hypothesis, is the cutoff value for the level of significance 
which as used in this analysis was taken as 0.05. 
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VI. ANALYSIS OP DATA 



The analog data v;ere digitized at 0.2 second intervals. 
The relationship betv/een the upper limit of the bandv;idth 
to be investigated, and the sampling interval At, is 



Af 

c 



1 

2At * 



( 6 . 1 ) 



Here f is the Nyquist frequency v/hich in this experiment 
c 

was established at 2.5 Hz. This was sufficiently high to 
avoid aliasing or folding problems. The maximum lag time 
was chosen as consistent at 10 percent of the record length. 
This results in twenty degrees of freedom for each computed 
spectra. The 80 percent confidence limits for tv.^enty 
degrees of freedom are between 0.71 and 1.60 of the measured 
spectral estimates according to the chi-squared distribution. 
Because the record lengths were different for each run, 
there was a different frequency resolution for each run. 

The frequency resolution for each run is shovm on the 
respective figures. 

A. VERTICAL VELOCITY MEASUREMENTS 

The predicted and measured vertical velocity energy- 
density spectra, and the phase angle, and coherence-spectra 
of the wave height and vertical velocity were computed for 
runs two through six at the depths shown in Table I. 
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1 . Predicted and Measured Vertical Velocity 
Energy-Density Spectra 



The vertical velocities were predicted using (4.^1). 
The predicted vertical velocity-spectrum v;as determined from 

The term In brackets Is the square of the transfer function, 
(4.6), relating the vertical velocity to p. 

Measured vertical velocity energy-density spectra, 
for Runs Tv;o and Six, are shovm In Figures (11) and (12). 

The principle energy band was from 0.04 to 0.37 Hz as deter- 
mined from the v/ave spectra. For Run Six, depth 5-70 meters 
the theoretical values of the vertical velocity spectrum 
were within three to five percent of the measured vertical 
velocity spectrum values over a frequency range of 0.04 
to 0.34 Hz. The true error of the velocities Is the square- 
root of the spectral differences and was approximately two 
percent. This is within the combined error of the v/ave and 
velocity meters; hence most of the measured energy-density 
within this bandwidth was v;ave Induced v;ith relatively 
little spectral energy-density being due to turbulence. 

On the other hand, for Run Number Two (the run of 
deepest depth at 14.62 meters) the absolute difference 
betv/een the measured and theoretical spectra was small over 
the much narrower bandwidth of 0.05 to 0.19 Hz. The large 
absolute difference between the measured and theoretical 
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Figure 11. Measured and Theoretical Vertical Velocity 
Spectra for Run Number Tv/o. Depth 1^1.62 Meters. 
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spectra at high frequencies for both runs indicates that 
beyond 0.19 Hz for Run Number Two and 0.3^ Hz for Run 
Number Six, that most of the energy was due to turbulence. 

The decrease of wave induced motion v;lth increased 
depth can also be seen by comparing the vertical velocity 
spectra for the two runs. For Run Number Six there v;as a 
peak in energy due to a three second wave component. This 
peak was not present in the spectra for Run Number Tvjo. 

That the energy present at this frequency and depth was 
due to turbulence can be seen by the large deviation 
betv;een the measured and theoretical spectra. 

Like the spectral energy-density in the wave Induced 
part of the spectrum, the turbulent energy-density level 
was one order of magnitude less for Run Number Tvjo compared 
to Run Number Six. This Indicates that turbulent motion 
decreased at the same rate as v;ave induced motion under the 
particular conditions of moderate sv.’ell, low winds, and 
intermediate water depth and that the Intensity of the 
turbulence is associated with the rms value of the 
viave-lnduced motion. 

A measure of the total energy available in the 
measured spectra is given by the variance. The value of 
the variance, and therefore the total kinetic wave-induced 
energy present, decreased v.'ith Increasing depth. There is 
almost an order of magnitude difference between the two 
values indicating a rapid decrease in energy vjith depth 
as predicted by theory. 



2 . Coherence of the VJave Height and Measured 

Vertical Velocity 

The coherence spectra are also shovm in Figures 
(11) and (12). In general ^ the coherence v;as 0.8 or 
greater throughout the range of significant vmve induced 
energy-density. This means that there is a high correla- 
tion between the vjaves and v/ave-lnduced vertical motion and 
that relatively little of the energy in the significant 
energy-density band was due to turbulence, but Instead, 
v;as v/ave-induced . 

3 . Phase Angle Between V/ave Height and the 

Vertical Velocity 

The phase angle theoretically leads wave height by 
90 degrees. The two computed phase spectra in Figures (11) 
and (12) further reflect the accuracy of linear theory at 
both shallow and deep depths. 

Bordy (1972) encountered the greatest departure from 
linear vjave theory in the computation of the phase angle. 
Much effort v;as expended during the course of this experi- 
ment to ensure that phase angle errors due to experimental 
technique were minimized. The concurrence of the constant 
phase angle vjith depth further supports linear theory. It 
is felt that keeping the wave sensor directly over the flow- 
meter was the sole reason for the difference in the results 
obtained in this research as compared with those of Bordy. 

. Comparison of Phase Angle and Coherence for 

Different Depths 

As can be seen from Figure (13) > there is a definite 



relation between the phase angle for a given frequency and 
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the depth of measurement. The 90 degree phase shift was 
observed to a period of 0.^ Hz at the shallov;est depth, 
but only as far as 0.25 Hz for the deepest case. 

The upper frequency limit of high coherence also 
decreases vjith increasing depth as Illustrated in Figure 
(l4). Both of these phenomena can be explained by linear 
theory. Linear theory predicts decreasing particle motion 
v/ith increasing depth due to a passing v/ave and that for a 
given depth the lov;er frequency vjaves induce greater parti- 
cle motion. This is reflected in the coherence and phase 
plots. The highest frequency at which the coherence 
approached one and at v;hich the phase compared vjell with 
theory for all analyzed runs corresponded v^ith the onset 
of significant turbulent energy contribution. 

B. HORIZONTAL VELOCITY MEASUREMENTS 

The predicted and measured total horizontal velocity 
energy-density spectra were computed for Run Number One. 

The phase angle and coherence spectra v;ere computed using 
the wave height and v-component of velocity. 

1. Predicted and Measured Horizontal Velocity Spectra 
The total theoretical horizontal velocity was 
computed from (^.3)* The theoretical energy-density 
spectrum of the total v;ave-induced horizontal velocity was 
calculated from the expression 



|u(f) 




cosh k(h+z)-,2 |„/^n| 2 
sinh kh ^ I 



(6.3) 



^9 




I 
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FREQUENCY (HZ) 

Figure 14. Vertical Coherence Spectra for Runs Two 
Through Six. 
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The term in brackets is (^.5)j the transfer function 
relating the total horizontal velocity to n • The total 
spectrum was found by calculating the individual u- and 
v-component spectra and summing over the frequency band 
of interest. The spectra are illustrated in semi-log format 
in Figure (15). 

The difference betv/een the spectra at the peak was 
found to be about 12 percent, giving an absolute value of 
the difference of about 3-5 percent in the significant 
energy-density frequency band. 

2 . Coherence of V/ave Heights and Measured Total 

Horizontal Velocity 

The coherence (of p and the v-component of velocity) 
was calculated to be greater than 0.9, Figure (15), through- 
out the significant energy frequency band of 0.03 to 0.07 
Hz. This high coherence was due to the fact that the v- 
component of velocity was primarily v/ave-induced and the 
turbulent contribution to the spectrum vjas small. 

3 . Phase Angle Betvjeen V/ave Height and the Total 

Horizontal Velocity 

The theoretical phase angle for the wave height 
and the v-component of velocity is zero degrees. This was 
found to be true in this case as illustrated again in 
Figure (15). 

C. PROBABILITY DENSITY FUNCTIONS AND STATISTICS OF THE 

DATA SETS. 

A probability density function v^as computed for each 
measured quantity for the six runs. The points comprising 
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Figure 15. Measured and Theoretical Horizontal Velocity 
Spectra for Run Number One. Depth 7.0A Meters. 
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a histogram of the data sets were plotted against both a 
Gaussian and a Gram-Charller distribution. The chl-squared 
goodness-of-f It test vras computed for both theoretical 
distributions. The statistics, variance, standard deviation, 
skevmess, and kurtosls were also computed for all eighteen 
data sets. These results are listed In Table IV. 

1 . Probability Density Functions 

A typical PDF Is shovm In Figure (l6). It Is 
knovm that the ocean's surface has approximately a Gaussian 
distribution. However, the non-llnearltles of the surface 
Introduce values of the skevmess and kurtosls that deviate 
from Gaussian and result In a Gram-Charller distribution. 

This Is explained by the physical fact that waves tend to 
form short, sharp crests and long troughs. For the waves 
this produces a positive skewness. Indicating more values 
In the time series belovj mean v;ater level than above. 

It follov/s that the v;ave Induced particle motion 
should follow the same behavior. Hovjever, of the eighteen 
data sets analyzed, three did not have positive skewness. 
These three were the velocity components u and w from Record 
Five and w from Record Three. Why this Is so, despite the 
fact that all the v;ave records were positively skev/ed, 
remains unexplained. A possible explanation could be that 
this peculiarity Is caused by Internal waves vjhlch were 
prevalent during the experiment. 
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TABLE IV 

CHI-SQUARED 

RUN/ STANDARD PARAMETER 

DEPTH (m) QUANTITY VARIANCE DEVISTION SKEWNESS KURTOSIS GAUSSIAN GRAM-CHAR 
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Figure 6. Gaussian and Gram-Charller Distribution Plotted Against Data PDF for 
the w-Component of Velocity Run Two. Depth 14.62 Meters, 8 June 1972. 



2. Goodness-of-Fit Test 



The chi-squared goodness-of-fit test was used to 
compare the computed PDF to both the Gaussian and Gram- 
Charlier distributions for all eighteen data sets. 
Qualitatively, based on the smaller chi-squared values 
obtained from a comparison of the PDF with the Gram- 
Charlier distribution, it can be stated that all eighteen 
records more closely fit the Gram-Charlier distribution 
than the Gaussian. The chi-squared parameter for all 
eighteen data sets is also listed in Table IV. 
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VII. CONCLUSIONS 



The measurements in this research vjere made under condi- 
tions of moderate swell and low wind velocity. V/ave height 
vjas measured vjith a penetrating v;ave staff and water 
particle motion measured with an electromagnetic flowmeter. 

A. ENERGY-DENSITY SPECTRA 

The wave-induced vertical velocity energy-density 
spectra were computed for six depths. The frequency of the 
peak energy-density remained constant vjlth depth. Hov;ever, 
the significant energy-density of the wave-induced motion at 
high frequency decreased with increasing depth according to- 
the attenuation aspects of linear theory. The magnitude of 
the energy-density spectra also decreased ^^^ith increasing 
depth according to theory. It was noted that the lowest 
frequency at which turbulence contributed significantly to 
the spectra shifted to lower frequencies v.’ith Increased 
depth. The intensity of the turbulent contribution decreased 
with depth and appears associated with the variance of the 
wave-induced velocities. 

B. COHERENCE 

Coherence of the wave height and wave-induced velocities 
was computed to be greater than 0.85 throughout the range 
of significant energy-density. The upper frequency limit of 
high coherence decreased with increasing depth. The value 
of the upper limit was established at the frequency where 



57 



turbulence became an Important factor in the energy-density. 
The near unity values of coherence over the band of signifi- 
cant spectral-density imply a near linear process betv;een 
vmve height and Induced vjater particle motion. It is felt 
that the small irregularities in coherence in the significant 
energy-density band were due to turbulence. 

C. PHASE ANGLE 

The computed phase angle spectra of the v;ave height and 
induced water particle motions were consistent with the 
theory at all depths for both vertical and total horizontal 
velocity cases. This is the part of the analysis v;here 
Bordy (1972) encountered the greatest deviation from theory. 
It is imperative that the v;ave staff be directly over the 
flov/meter in order that there be no phase error. A similar 
decreasing upper frequency limit with increasing depth for 
consistent phase angle values v;as present as it v;as for both 
the coherence and the energy-density spectra. 

D. STATISTICAL INTERFERENCES 

The probability density functions of all measured 
quantities were qualitatively found to be better described 
by the Gram-Charlier distribution than by the Gaussian using 
the chi-squared goodness-of-fit test. The variance of the 
vertical velocities, a measure of the kinetic energy present, 
decreased v/ith increasing depth vjhich is consistent with 
linear theory. The skewness of the wave height frequency 
distribution v;as found to be positive in all cases as a 
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consequence of the fact that vmves tend to have short, 
sharp crests and elongated troughs. The values of skewness 
with fev; exceptions, v;ere found to be positive for the v;ave 
induced quantities as v;ell. These exceptions might be 
explained by the presence of internal waves. 
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